Cerenkov radio pulses from electromagnetic showers in the time-domain 
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The electric field of the Cerenkov radio pulse produced by a single charged particle track in a 
dielectric medium is derived from first principles. An algorithm is developed to obtain the pulse 
in the time domain for numerical calculations. The algorithm is implemented in a Monte Carlo 
simulation of electromagnetic showers in dense media (specifically designed for coherent radio emis- 
sion applications) as might be induced by interactions of ultra-high energy neutrinos. The coherent 
Cerenkov radio emission produced by such showers is obtained simultaneously both in the time and 
frequency domains. A consistency check performed by Fourier-transforming the pulse in time and 
comparing it to the frequency spectrum obtained directly in the simulations yields, as expected, 
fully consistent results. The reversal of the time structure inside the Cerenkov cone and the signs 
of the corresponding pulses are addressed in detail. The results, besides testing algorithms used for 
reference calculations in the frequency domain, shed new light into the properties of the radio pulse 
in the time domain. The shape of the pulse in the time domain is directly related to the depth 
development of the excess charge in the shower and its width to the observation angle with respect 
to the Cerenkov direction. This information can be of great practical importance for interpreting 
actual data. 
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I. INTRODUCTION 

It was nearly 50 years ago that Askaryan proposed 
to detect high energy particles through the coherent 
pulse they emit as they interact in a dense medium 
As secondary electrons, positrons and gamma rays are 
produced they develop electromagnetic showers in the 
medium which acquire an excess negative charge, which 
Askaryan estimated to be of order 10% of the total num- 
ber of electrons and positrons. This is so in spite of the 
interactions being completely charge symmetric, because 
matter in the medium only contains electrons. M0ller, 
Bhabha and Compton scattering of matter electrons, ac- 
celerate them into the shower while electron-positron an- 
nihilation and Bhabha scattering decelerate the shower 
positrons thus also contributing to the excess charge, a 
mechanism referred to as the Askaryan effect. A more ac- 
curate calculation of the Askaryan effect indicated that 
the excess charge is actually ~ 25% of the total number of 
electrons and positrons [1] . Such an excess charge devel- 
ops a coherent electromagnetic pulse as it travels through 
a non absorptive dielectric medium. The coherent part 
of the pulse is mainly due to the wavelength components 
which are large compared to the shower width. The en- 
ergy radiated in the coherent pulse scales with the square 
of the excess charge and hence with the square of the 
shower energy. Such scaling naturally makes the detec- 



tion of coherent radio pulses an attractive and promising 
technique for the detection of ultra high energy particles, 
such as cosmic rays. 

Radio detection of air showers was extensively studied 
in the 60's and 70's f^. The drive to detect high energy 
neutrinos in the late 80's turned back the attention onto 
radio pulsesproduced by them in dense media such as 
natural ice [3] or the regolith beneath the Moon's sur- 
face [1] . The first full simulations of the Askaryan effect 
and the coherent pulses created in dense media were ob- 
tained in the ear ly 90's which allowed more quanti- 
tative calculations and experimental programs were soon 
after started to search for neutrinos with arrays of an- 
tennas at Antarctica Q or with radio telescopes from 
Earth \S\ . The Askaryan effect was measured for the first 
time firing photon bunches into sand at SLAC in 2000 M 

- and later in other dielectric media including ice 

- and since then the field has received an enormous boost, 
strengthening previous initiatives using antennas buried 
in ice [H, and radiotelescopes , and developing 
new ones such as a balloon flown antenna array |16l4l8l |. 
new radiotelescope searches [T9l - [2^ and new radio mea- 
surements of air showers 2^. 

The first calculation of the radio emission from elec- 
tromagnetic showers used a specifically designed Monte 
Carlo simulation code - the ZHS code - to calculate co- 
herent radio pulses in ice 0, Q . The code has been ex- 
tended to include the LPM effect [231, to calculate in an 
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approximate manner hadronic showers [26j and neutrino- 
induced showers [l^l , to treat other dielectric media [11] , 
and to perform an optimal statistical thinning that al- 
lows the simulation of pulses from ultrahigh energy show- 
ers [2^, and remains as a reference in the field. This 
code was designed to calculate the Fourier components 
of the electric field in the frequency domain. Alternative 
simulations using other codes such as GEANT3 [30, 31], 
GEANT4 [MQ and the AIRES-HTIERRAS [H HI 
code, have yielded results compatible to within ^ 5%. 
Semi-analytical calculations have also been performed 
[36j . All of these use the same technique to calculate the 
radio pulse in the frequency domain, but to our knowl- 
edge no full calculation exists in the time domain yet. 

All experimental arrangements measure the electric 
field as a function of time, and full understanding of the 
properties of the pulse as a function of time is thus also 
very important. Although the conversion from the fre- 
quency to the time domain is in principle straightforward 
and the algorithm in ZHS computes all required informa- 
tion to obtain it, there have been a number of doubts con- 
cerning the unconventional choice of Fourier transform as 
used in the code 0] , as well as the sign, phase and causal- 
ity properties of the pulse [1^ , that have complicated the 
analysis and interpretation of data. 

In this article we develop a formalism to calculate the 
pulse directly in the time domain. We simultaneously 
calculate the pulse of the same electromagnetic shower 
in both the time and frequency domains. An exhaustive 
comparison yields fully compatible results, makes patent 
the relative advantages of each approach, and sheds new 
light into the properties of the radio pulse in the time 
domain which can be related to those of the shower and 
can be of great practical importance in interpreting ac- 
tual data. Some of these properties are discussed in more 
detail suggesting possible applications. 

Although the method developed in [^l, and extended 
here to the time-domain, has been obtained in the frame- 
work of Cerenkov radiation, it derives directly from 
Maxwell's equations and addresses classical radiation 
from charges in a pretty general fashion. Simple ex- 
tensions of this work can be used for instance to calcu- 
late transition radiaton as particles cross dielectric me- 
dia interfaces or to calculate the complete radiation pat- 
terns from charges moving in magnetic fields including 
Cerenkov radiation, as has been known for long to be 
important for ultra high energy air showers. 

This paper is structured as follows. In Section [II] we 
rederive the expression for the electric field in both the 
time and frequency domain in a form that can be easily 
used for practical applications and make the connection 
to the expression derived in the original ZHS paper 
We also discuss some simple current density models and 
relate them to the results of a full electromagnetic shower 
simulation. In Section iHll we perform a consistency check 
by Fourier-transforming the pulse in time and comparing 
it to the frequency spectrum obtained in the simulations. 
The summary and outlook constitute the last section. 



II. THEORY AND MONTE CARLO 
IMPLEMENTATION 

A. Theory 

We start from Maxwell's equations for linear, isotropic, 
homogeneous and non dispersive media. In the Interna- 
tional System of units: 



V • E : 



V -B = 



V X E = -■ 



dB 

dt 



(1) 



aE 



V X B = + ne— (2) 
ot 



where p is the charge density of the source, e = ei Co and 
fi = /Zr/xo are the total permittivity and permeability ex- 
pressed in terms of the relative (/x^ and e^) and free space 
(/io and eo) permittivities and permeabilities. All effects 
of induced currents and electric polarization are automat- 
ically accounted for by the displacement field D = eE 
proportional to the electric field, E and the magnetic 
field strength H = {fj,)~'-^'B, proportional to the magnetic 
field, B. 

We recall the formal solution introducing the vector 
and scalar potentials (A and (p): 



B = V X A 
OA 



E 



(3) 
(4) 



that naturally satisfy V • B = 0, and the equation involv- 
ing the V X E term. Choosing the transverse gauge, in 
which V • A = 0, the two remaining equations imply: 



V^A 



6>2A 



P 

€ 



(5) 
(6) 



where Jj^ is the transverse current, a divergenceless com- 
ponent of the current density, which in the limit of obser- 
vation at large distances from the source can be shown to 
correspond to the projection of the current density per- 
pendicular to the direction of observation (of unit vector 
u), i.e., J_L = — u X (u X J). 

Both equations can be formally solved using Green's 
functions: 



1 
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|x-x'| 



(7) 



J±(x',t 



^(5 (V^elx - x'l - (t - t')) d^^dt' 



(8) 



The first is the familiar solution from electrostatics for 
the potential produced at the position x by a source with 
charge density p{x' ,t'). The second is the solution of the 
wave equation with wave velocity (eo/J.oer/J.r)~^ smaller 
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than the velocity of light in vacuum, c = (eo/^o) ^ i by a 
factor n = (erAir)~2 , the index of refraction. The Green's 
function for the wave equation involves a delta function 
that gives the familiar retarded time, t\ earlier than the 
observing time t. To evaluate the field at time t at a 
given position x, the current is to be evaluated at a time 
retarded by the time taken by light to reach observation 
point from point x', i.e. |x — x'|n/c. 



In this equation we can take a = (1 — npcosO) and 
X — t—nR/c—{l — n(3cos9)ti 2 which allows us to rewrite 
Eq.® as: 



A=J^ 
AttR 



(1 - nPcose)ti) - e(t - ^ - (1 - n/3 cos 61)^2) 



{l-n(3cose) 



(12) 



Radiation from charges traveling in straight 
lines 



We consider the shower as a superposition of finite par- 
ticle tracks of constant velocity. Each track is completely 
defined by two limiting times ti and its velocity v and 
the position vector of an arbitrary point of the track, xq , 
which we have chosen to correspond to the time t = 0. 
The transverse current density entering in Eq. ([5]) for a 
point charge moving with constant velocity, v, between 
the two end points simply reads : 



J_L(x',t') = ev^^^ (^yJ 



xo - vt') [Q{t' - h) - e{t' - t2)] 

(9) 

where — e is the charge of an electron, is the projection 
of the velocity onto a plane perpendicular to the direction 
of observation (recall that we consider large distances so 
that this direction is uniquely defined), and Q{x) is the 
Heaviside step function. 

We can now substitute the transverse current into 
Eq.®, integrate the three dimensional delta function 
substituting x' for xq + vt' and approximate the dis- 
tance between x and x' by |x — Xq — vt'\ ~ i? — v • ut' , 
where we define i? = |x — xo|. In the limit of large dis- 
tances of observation the denominator |x — x'| can be 
simply approximated by R. However, we must use the 
above approximation in the argument of the retarding 
delta function to account for interference effects. This 
corresponds to the Fraunhofer approximation, in which 
the path difference between light pulses emitted at points 
Xq and x' = Xq + 'vt' is simply the distance between them 
projected onto the direction of observation. As a result 
the delta function reads S (t'{l — n/3cos0) — (t — 
with V — f3c, which can be cast into: 



1 



\l-n/3cose\ 



5 t' 



c 

1 — n/3 cos i 



(10) 



We note that the recurring factor (1 — nf3 cos 9), with 9 
the angle between v and u, gives zero for the Cerenkov 
angle 9c ■ Moreover the factor changes sign from positive 
to negative as the observation angle changes from being 
larger to smaller than the Cerenkov angle. Now we can 
perform the integration in t' in Eq. ([8]) which simply im- 

plies replacing t' in the step functions by ]^_„^ cosfl - 
now make use of the fact that: 



e{ax) 



Q{x) if a > 0, 
l-e(a;) ifa<0 



(11) 



Note that the modulus in the denominator of Eq. (ITOl) is 
removed because of an effective sgn(l — n/3 cos 0) that 
appears when changing the argument in the two step 
functions (according to Eq. (|lip ). This expression is easy 
to implement in a shower simulation by splitting particle 
tracks in portions that can be approximated by uniform 
motion. 

As 9 approaches the Cerenkov angle 9c the numerator 
and denominator of Eq. (ll2|) approach zero. To obtain 
a formal limit for the Cerenkov angle we multiply and 
divide by St to obtain: 



i?A(t, 



47reoC' 



e {t ~ ^ - {1 - 71/3 cos 9)ti) ~Q{t - ^ - {l-nl3cos 9)t2) 



(1 - nl3cos9)5t 



(13) 



The limit — >■ 6'(7 is equivalent to (1 — n(3 cos 9)6t — > 
which can be shown to give the first derivative of the step 
function, the function S{t). The limit is then: 



RA{t,9c) 



eflr 



6{t-'^)v^St (14) 



We note that the vector potential in this limit is simply 
proportional (and parallel) to ^±61, i.e. to the projection 
of the displacement vector onto a plane perpendicular to 
the observation direction. This expression can also be 
implemented in a shower simulation provided a suitable 
approximation is made for the delta function. 

The expression for the electric field is given by Eq. (|4]) 
and only the term with the time derivative of the vector 
potential gives contribution to the radiation term so that: 



R'E{t,9) = - 
S{t 



^ - {l-nl3cos9)ti 



nR 



(1 — nf3 COS 9)t2 



(1 - n/3 COS 6*) 



(15) 



The radiation field due to a single particle track with 
e > is similar to the schematic diagram shown in Fig. [T] 
Such a particle produces radiation when the track starts 
or ends. The two pulses "as seen" by the observer (placed 
at angle 9 w.r.t. the particle track) are separated by a 
time interval associated to the difference in propagation 
time (1 — nP COS 9)St. 

Let us first consider an angle exceeding the Cerenkov 
angle so that (1 — n/3 cos 9) is positive. The electric field of 
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FIG. 1. Schematic representation of the radiation fields in the 
time domain induced by a single particle with positive charge 
e > traveling at constant velocity (5 between times ti and 
t2. Top panel vector potential (see Ea. p2p ). Bottom panel 
electric field (see Ea. (|15p ). See also text for more details. 



the first pulse corresponds to the start point of the track 
{ti) and it is anti-parallel to vj. according to Eg. ([T5l) . 
while it is parallel for the second pulse which corresponds 
to the end point (^2)- The sign of the electric field pulse is 
opposite to the sign of the particle acceleration in both 
cases. The zero in the shown arrival time is arbitrary 
and corresponds to t = nR/c, i.e. it is a reference time 
associated to the arrival of a signal from the reference 
position xq. The two pulses associated with the track 
take place later than this reference time. 

As the angle decreases and becomes smaller than the 
Cerenkov angle, the situation is reversed: The first pulse 
corresponds to the end point of the track (^2), while the 
second corresponds to the start point (^i ) . Moreover not 
only is the arrival of the pulses as seen by the observer 
inverted, but both take place before the reference time. 
This apparent acausal behavior is due to the fact that the 
particle travels at a speed greater than that of light in the 
medium. Although the terms responsible for the first and 
second pulses are interchanged, and there is a sign change 
associated with this interchange, it is compensated by 
the denominator of Eq.(|31) that also reverses its sign. As 
a result there is no change in the sign of the electric 
field of the first and second pulses as the Cerenkov angle 
is crossed, and the double peak structure at any given 
time has the same qualitative behavior as the observation 
angle changes. This seems physically sound since there 
can be no discontinuity of the electric field across the 



Cerenkov cone boundary. 

For observation at the Cerenkov angle both signals ar- 
rive simultaneously. In this limiting case the electric field 
can be formally obtained taking minus the derivative of 
the delta function given by Ea.(IT4|. This again corre- 
sponds to a double pulse first antiparallel and then par- 
allel to Vj_. 



C. Equations in the Frequency Domain 

The expression for the electric field in the frequency 
domain used in the ZHS simulation code (Eq.(12) in f3|) 
reads: 



E(w,x) = 



e/ir e 



ikR 



R 



gi(w-k-v)ti ^ 



— k- v)(5t Y 



i(LO — k • v) 
(16) 

We recall that this equation has been obtained with the 
following convention for the Fourier transform of the elec- 
tric field: 



/» = 2 / f{t) e^'^'dt 



(17) 



where the factor 2 corresponds to an unusual convention 
(this factor is usually either 1 or (27r)^2). Applying this 
Fourier transform definition to Eq. (1151) giving the electric 
field in the time domain we obtain: 



E(c^,x) 



1 

— V I 



27reoc2 R 

giLj[n-R/c+(l — n/3 cos 6)]ti 



piLu[nR/c-\-{l—n(3 cos d)]t2 



(1 — nl3 cost 



(18) 



which can be easily rearranged to give exactly Eq. p6)) 
noting that k = Moreover if we apply the Fourier 
transform to Eq.(fH|) which applies to the limit 9 Oc 
we get: 



i?A(a;,x) 



(19) 



The electric field is obtained taking minus the time 
derivative which in Fourier space is just a factor iu), giv- 
ing again the same result as Eq.(13) in for the electric 
field in the frequency domain at the Cerenkov angle. 

These calculations show the consistency of Eq. (ITSI) ob- 
tained in the time domain with Ea.(|16p which gives the 
field in the frequency domain: They are simply Fourier 
transforms of each other as expected. 



D. Pulses for Simple Charge Distributions 

Before performing a Monte Carlo simulation of electro- 
magnetic showers, it is interesting to extend the calcula- 
tions to simple models for the shower. These models al- 
low us to obtain relations between the shape of the pulse 
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in the time domain and the time and spatial distribution 
of the charge. 

A simple yet interesting model consists of a charge 
Q{z') that rises and falls along the shower direction z' 
and spreads laterally in x' and y' . Assuming cylindri- 
cal symmetry we can write the current associated to this 
charge distribution as: 



J(x',t') = vf{z',r')Q{z')5{z' - vt') 



(20) 



Here r' is a two dimensional vector in the {x',y') plane 
transverse to z\ and the function f{z',r') gives the 
charge distribution in such a plane as a function of shower 
depth, with normalization chosen so that Q indeed gives 
the excess charge: 



J dhfiz',r) 



2tt 



/(z',r',</>') = l (21) 



with (f)' the azimuthal angle in cylindrical coordinates. 

The simplest case is that of a line current along the z' 
direction without lateral extension in which f{x', y') is re- 
placed by the two dimensional delta function 5{x')5{y'). 
This approximation was also discussed in [27| in the fre- 
quency domain, where it was referred to as the one- 
dimensional approximation. When such line current is 
substituted into Eq.(l8]) and integrated in x', y' and t' 
making the Fraunhofer approximation, a relatively sim- 
ple expression is obtained that relates the vector poten- 
tial in the time domain to the excess charge Q{z'): 



dz'Q{z')5 



z'{l - nf] cose) - V (t 



nR 



(22) 



The delta function relates the depth in the shower de- 
velopment z' to the observation time t through a linear 
function: 



z' = at) = p 



ct — nR 
1 — n/3 cos ( 



(23) 



As the observation angle approaches the Cerenkov angle, 
the time interval corresponding to the depth spanned by 
the shower, i.e. the pulse width, becomes smaller. We 
thus recover a familiar result already discussed in al- 
though in the frequency domain. 

Performing the integration in Eq. ([22|) yields. 



RA 



Iicf3 



47r |1 - 71/3 cos 6\ 



Qm) 



(24) 



The factor {l — nf3 cos9)~^ arises from applying the chain 
rule to the derivative of Q[({t)]. As a result the pulse in 
the time domain can be regarded as the derivative of 
the development of the charge excess along the shower, 
scaled with the Cherenkov factors (1 — n/S cosd)^^ and 
|1 — n/3cos0|~^, and converted from depth into time 
through Eq. (P5)) . i.e., the pulse is firstly positive and then 
negative with respect to vj_ since in a real shower Q{z') 
corresponds to an excess of negative charge. 

A number of interesting results can be directly read off 
Eg. ipS)) . If the development curve for the excess charge 
Q{z') is not symmetric, as happens in real showers, the 
asymmetry in its derivative is directly reflected into an 
asymmetry between the negative and positive parts of 
the pulse. Also it is interesting to note that when the 
angle of observation is below the Cerenkov angle, the 
pulse shape is inverted in time because the early part 
of the pulse corresponds to the end of the shower while 
the beginning of the shower corresponds to the end part 
of the pulse, as explained above. Still, the polarity of 
the first and second pulses remains the same because, 
although the slopes before and after shower maximum 
change sign, there is an extra sign change induced by the 
factor (1 — n/3 cos 0)~^ . This is in complete analogy to 
what was discussed for a single track. 

In the case of observation in the Cerenkov direction, 
the z' dependence of the delta function in Eq.([22|) disap- 
pears and the delta function can be factored away from 
the integral, to give a pulse of amplitude directly propor- 
tional to the integrated excess track length of the shower. 
The delta function term is due to all parts of the line 
current being observed simultaneously at the Cerenkov 
angle. These are two familiar results already emphasized 
in [2]. 

The simulation has shown that the model with the ab- 
sence of a lateral distribution breaks down at \0 — 6c\ ^ 
2.5°. This result is consistent with that found in [2^ 
where the one dimensional model was studied in the fre- 
quency domain. 

It is instructive to extend the line current model to 
a more realistic three dimensional current /(z',r) with 
cylindrical symmetry and current given in Ea. (j20p . In 
that case the expression for the vector potential with 
two delta functions can be integrated in t' and (p' and 
the resulting expression involves a double integral over 
cylindrical coordinate r' and the shower depth z': 



RA = v_L 



2tt 



r'dr' 



dz'f(z',r')Q{z') 



Q{nPr' sine - \z'{l - n/3 cos 0) - {vt-nl3R)\) (26) 
^J[nPr' sinef - [z' {I - n^ cose) - {vt-n(3R)f 



where the delta function in Eg. ip^ introduces a factor 
\l-nPcose\-'^. 

The electric field is obtained taking minus the deriva- 
tive of the vector potential with respect to time: 



RE= - 



lic(3 



dQiO 



47r (l-n/3cos6')|l-n^cos6'| dC 



(25) 



This expression, despite being more cumbersome than 
Ea. (l22)) . if solved analytically for realistic lateral distribu- 
tion functions, could give insight into useful parametriza- 
tions of the pulse in the time domain. In any case it can 
be used for numerical simulations. 
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In the Cerenkov hmit Ea. ((26l) becomes 
^1 w±_ 



RA = 



2tt nP sin 9c J- 

r'dr'f(z',r') 



dz'Q{z') 



(27) 



vt-n/3R 
n/3 sin &c 



This equation shows that the non-zero width of the elec- 
tromagnetic pulse at the Cerenkov angle is the result of 
the lateral distribution of the shower. Although the inte- 
gral is rather complicated to evaluate for realistic lateral 
shower profiles, it can be shown that for distributions of 
the form f{r') = (r')~" for integers n > 2 the electric 
field E oc vxsgn(t — nR/c)\vt — n/3R\^", which is a fast 
bi-polar pulse of non-zero width. 

This model still has some limitations. Note that 
Ea.(|27p predicts a pulse that is symmetric in time while 
simulations have shown that the pulse at the Cerenkov 
angle is asymmetric. This is due in part to the radial dis- 
tribution of velocities of the shower which is not included 
in the model. The development of a current density vec- 
tor model that can accurately produce the features of 
Cerenkov radiation is work in progress. 



E. Implementation in the ZHS Monte Carlo 

The ZHS Monte Carlo [2| allows the simulation of 
electromagnetic showers and their associated coherent 
radio emission up to EeV energies [2£]. Originally de- 
veloped in ice Q, it has been extended so that electro- 
magnetic showers in other dielectric homogeneous me- 
dia can be simulated [H, [2^. The code accounts for 
bremsstrahlung, pair production, and the four inter- 
actions responsible for the development of the excess 
charge, namely M0ller, Bhabha, Compton scattering and 
electron-positron annihilation. In addition multiple elas- 
tic scattering (according to Moliere's theory) and contin- 
uous ionization losses are also implemented. The elec- 
tron/positron tracks between each interaction are split 
into subtracks so that no subtrack exceeds a maximum 
depth fixed at 0.1 radiation lengths. For low energy par- 
ticles these subdivisions are actually reduced to ensure 
that no subtrack is comparable to the particle range, and 
they become the step used to evaluate ionization losses 
and multiple elastic scattering. Convergence of results as 
the step is reduced has been carefully checked fssj . 

In order to account for interference effects between the 
radiation emitted due to the particles responsible for the 
excess negative charge, the ZHS code was designed to fol- 
low all electrons and positrons down to 100 keV kinetic 
energy threshold, as well as to carefully account for time, 
by considering deviations with respect to a plane front 
moving at the speed of light injected in phase with the 
primary particle. In addition to the delays associated to 
the propagation geometry, those due to particles trav- 
elling at velocities smaller than the velocity of light are 



accounted for assuming the energy loss is uniform across 
the step. An approximate account is also made of the 
time delay associated to the multiple elastic scattering 
processes along the step. 

As a result the tracks of all charged particles in a 
shower are divided into multiple subtracks which are as- 
sumed to be straight and to have constant velocity. The 
positions of the end points of these subtracks as well as 
the corresponding times are readily available by design, 
and they can be used to compute the frequency compo- 
nents of the electric field making extensive use of Eq. ([T5| , 
taking into account the relative phase shift between dif- 
ferent tracks because of their different starting point po- 
sitions and time delays. 

In this work we have extended the Monte Carlo to also 
calculate the pulse in the time domain. A routine has 
been developed to account for contributions of each of 
these particle subtracks to the vector potential, making 
extensive use of Ea.([T^. Each subtrack contributes a 
unit "rectangle" to the vector potential, which varies in 
height, "duration" and sign - see Fig. [U depending on 
the velocity, the relative orientation of the track with 
respect to the direction of observation and the charge of 
the particle. When the observation direction is very close 
to the Cerenkov angle the delta function in Eg. dTi]) is 
replaced by a rectangle corresponding to a nascent delta 
function j37j . If the sampling time bin width is set to AT 
then a natural choice of nascent delta function is given 
by 



1 AT ^ J. < AT 

AT' 2 ^ - 2 

0, otherwise 



(28) 



In this case the base of the rectangle is fixed by the intrin- 
sic "time resolution" AT of the simulation and the pulse 
height depends on AT. In practice, the time domain ra- 
dio signal can be reconstructed with an antenna receiver 
system and digital sampling electronics. The time reso- 
lution of a single waveform is determined by the digital 
sampling bin width and the high frequency cutoff of the 
receiver system. 

Once the vector potential induced by each subtrack is 
defined, the contribution of all charged subtracks in the 
shower is obtained and the vector potential is derived 
with respect to time to obtain the electric field in the 
time domain. 

In the next Section we show several examples of the 
results of this procedure. 



III. RESULTS 

In Fig. [5] we show the electric field as a function of 
arrival time of the signal obtained with the ZHS code in 
a single 1 PeV electron-induced shower in ice for different 
observation angles. The zero in the shown arrival time 
is measured with respect to the arrival time of a pulse 
emitted as the primary particle initiating the shower is 
injected in the medium. 
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FIG. 2. Electric field as a function of time as obtained in 
ZHS simulations of a single 1 PeV electron-induced shower in 
ice for different observation angles. Top panel: Observation 
at the Cerenkov angle, bottom panel; observation at 9c — 5° 
(long green dashes) and at 9c + 5° (short blue dashes). In 
the bottom panel the red solid histograms represent the elec- 
tric field obtained applying Ea. (|25p to the simulated excess 
negative charge Q(z). 



FIG. 3. Top panel: Longitudinal development of the excess 
negative charge as obtained in ZHS simulations of 1 PeV 
(long green dashes) and 100 PeV (short blue dashes) electron- 
induced showers in ice. Bottom panel: Electric field as a func- 
tion of time generated in the showers shown in the top panel 
(dashed histograms), for observation angle 9c + 10°. The 
solid histograms represent the electric field obtained applying 
Ea. (|25|) to the simulated excess charge Q{z) shown in the top 
panel. 



The electric field is parallel to the projection of the 
velocity onto a plane perpendicular to the direction of 
the observation at early times and anti-parallel later on. 
This is expected after the discussion in Section lll Bl of the 
electric field emitted by a single positively charged par- 
ticle, with the important difference that in a shower the 
electric field is produced by an excess of negative charge 
and the polarity of the field is reversed with respect to 
that shown in Fig. [TJ Also as in the case of a single track 
there is no change in the polarity of the pulse when ob- 
serving inside {9 < 9c) or outside {9 > 9c) the Cerenkov 
cone {9c)- The pulse always starts being positive (paral- 
lel to Vj_) and ends being negative (antiparallel to vj^) 
regardless of the observation angle. This feature can be 
used as a discriminator against background events for 
neutrino searches. It can be also clearly seen that the 



pulse is broader in time away from the Cerenkov cone 
than close to it with an apparent duration proportional 
to Az\l — n(3 cos 9\/c with Az being the spread along the 
shower axis of the excess charge (see Eg. ([25)1 V For ob- 
servation at the Cerenkov angle the apparent duration of 
the pulse is not zero, despite the fact that the Cerenkov 
factor 1 1 — nf3 cos 9c\ — > 0, because the shower spreads out 
also in the lateral dimensions {x and y directions). Also 
due to our definition of t = and to the presence of the 
Cerenkov factor in the (5— functions in Eq.([T5]), the pulse 
occurs at t > outside the Cerenkov cone and at t < 
inside it. 

According to the simple model developed in Section 
IIIDI the field away from the Cerenkov angle is propor- 
tional to the derivative of the excess charge distribution 



Q{z) with respect to t - Ea. ([25|) - or equivalently the 
derivative with respect to z since there is a Unear relation 
between t and z - Eg. (f23|) . The ZHS code also gives the 
longitudinal profile of the excess charge and we have ap- 
plied Eg.pSj) to the simulated Q{z), and compared to the 
electric field obtained directly in the Monte Carlo. This 
is also shown in Fig. [51 The agreement between the elec- 
tric field obtained directly in the Monte Carlo simulation 
(dashed histograms) and what is predicted by Eg. ([25)) 
(solid histograms) is remarkable. The electric field fol- 
lows the variation of the excess charge in z or equivalently 
in t. This explains why for a fixed observation angle the 
pulse changes sign from early to late times (for a typical 
shower Q{z) grows relatively fast, reaches a maximum, 
and then decreases more slowly with depth), and why 
it is asymmetric with respect to the time axis {Q{z) is 
not a symmetric function around its maximum). Also 
when the direction of observation is inside the Cerenkov 
cone, the observer sees the derivative of the beginning of 
the excess charge distribution first and the corresponding 
derivative of the end oi Q{z) at later times, while the op- 
posite is true for observations outside the Cerenkov cone. 
As a consequence the pulse at 6 < 9c looks like an an- 
tisymmetric copy with respect to t = of the pulse at 
6 > 9c, as can be clearly seen in Fig. [21 An accurate 
reconstruction of the time domain electric field could in 
principle determine on which side of the Cerenkov cone 
the event was observed. On the other hand the shape 
of the pulse can be conversely used to infer the depth 
development of the shower. 

Ea. ([^ stresses the fact that the features of the excess 
charge distribution are "mapped" in the time structure 
of the pulse. In particular it is well known that electro- 
magnetic showers with energies above the energy scale 
at which the LPM effect [331 starts to be effective (~ 
PeV in ice [111), are "stretched" in the longitudinal di- 
mension and often show peaks in their profile [251 |40| - 
l42j . These two features should translate into the dura- 
tion in time also of the pulse and into its time structure 
that should also exhibit multiple peaks. This is shown in 
Fig. [21 in which due to the LPM effect the longitudinal 
profile of a 100 PeV electron-induced shower exhibits two 
peaks which appear as 2 positive and 2 negative peaks 
in the time structure of the pulse. For comparison a 1 
PeV electron-induced shower not affected by the LPM 
effect and its corresponding electric field are also shown. 
The linear relation between the time domain structure of 
an electric field and the shower profile suggests that the 
longitudinal profile of the shower could be reconstructed 
from an observation off the Cerenkov angle. 

The extended ZHS code is able to calculate both the 
electric field as a function of time and its Fourier trans- 
form from first principles. Moreover, the two calculations 
can be made simultaneously for the same shower. Both 
calculations can be easily compared by performing the 
Fourier transform of the pulse calculated in the time do- 
main, following the convention in Ea. (jl7p . This provides 
a further check of the two methods, as well as a test of ac- 




Frequency [MHz] 

FIG. 4. Electric field frequency spectrum obtained in ZHS 
simulations of a single 1 PeV electron- induced shower in ice for 
different observation angles (green dashed lines). Also shown 
is the Fast Fourier Transform (FFT) of the electric field in the 
time domain obtained in simultaneous ZHS simulations of the 
same shower for two different time resolutions AT = 0.1 ns 
(red solid lines) and AT — 0.5 ns (magenta dotted line - only 
shown at the Cerenkov angle for clarity). 

curacy in the numerical procedures involved in the calcu- 
lation of the radio emission in both domains. An example 
is shown in Fig. [31 where the electric field as a function 
of frequency as obtained in ZHS simulations of a single 
1 PeV electron-induced shower is plotted along with the 
(Fast) Fourier Transform (FFT) of the electric field in 
the time domain obtained in simultaneous ZHS simula- 
tions of the same shower. The agreement between both 
spectra is very good for frequencies below wat ^ 2n/AT 
with AT an arbitrary time resolution needed for the ZHS 
simulations in the time domain. We do not expect to be 
able to reproduce the frequency spectrum at frequencies 
above ojat - proportional to the Nyquist frequency of the 
system. To illustrate this point in Fig. [H we also show 
the Fourier transformed spectrum (at the Cerenkov an- 
gle) of several time domain calculations performed with 
different time resolutions AT — 0.1 and 0.5 ns. One can 
see that the agreement between the frequency spectrum 
obtained in ZHS and the Fourier transformed time do- 
main electric field improves as AT decreases as expected. 
Calculations in the frequency domain are more advisable 
near the Cerenkov angle. 
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IV. SUMMARY AND OUTLOOK 

In this work we have developed an algorithm to obtain 
the Ccrcnkov radio pulse produced by a single charged 
particle track in a dielectric medium. We have imple- 
mented this algorithm in the ZHS Monte Carlo with 
which we can predict the Cercnkov coherent radio emis- 
sion emission of electromagnetic showers in dense dielec- 
tric media in both the time and frequency domains. 

An observer in the Fraunhofer region, far from the axis 
of the electromagnetic shower at an angle 6, sees a bi- 
polar pulse due to the excess of negative charge in the 
shower. The apparent time duration of the pulse is pro- 
portional to A2: (1 — n cos 9)/c with Az the spread of the 
shower in the longitudinal direction. At the Cercnkov 
angle (1 — n cos 6c) — > and the duration of the pulse 
is mainly determined by the lateral extent of the shower. 
At angles 6 > 6c, the observer sees first the electric field 
produced by the early stages of the shower, and the field 
due to the end of the shower later on, while the time se- 
quence reverses for observation at ^ < ^c- Regardless of 
the observation angle, the bulk of the electric field due 
to the excess negative charge is directed along v^ - the 
projection of the particle velocity onto a plane perpen- 
dicular to shower axis - at early times and in the opposite 
direction later on. The shape of the pulse maps the vari- 
ation with depth of the excess charge in the shower. This 
information can be of great practical importance for in- 
terpreting actual data. 

A consistency check performed by Fourier- 



transforming the pulse in time and comparing it to 
the frequency spectrum obtained directly in the simu- 
lations yields, as expected, fully consistent results. Our 
results, besides testing algorithms used for reference 
calculations in the frequency domain, shed new light 
into the properties of the radio pulse in the time domain. 

In the future we plan to implement the algorithm 
for time-domain calculations of electric field pulses in 
Monte Carlo simulations of hadronic and neutrino- 
induced showers, of great importance for neutrino de- 
tectors using the radio Cerenkov technique. Also we will 
explore how actual experiments can exploit the richness 
of information contained in the shape in time of the radio 
pulse to obtain information on the shower development. 
This could be of great help in the reconstruction of the 
parameters of the neutrino- induced showers and to dis- 
criminate against background events. 
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